Josephson current transport through a Quantum Dot in an Aharonov-Bohm Ring 
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The Josephson current through an Aharonov-Bohm (AB) interferometer, in which a quantum dot 
(QD) is situated on one arm and a magnetic flux $ threads through the ring, has been investigated. 
With the existence of the magnetic flux, the relation of the Josephson current and the superconductor 
phase is complex, and the system can be adjusted to ir junction by either modulating the magnetic 
flux or the QD's energy level Ed- Due to the electron-hole symmetry, the Josephson current I has 
the property I(ed,&) = /(— £d, $ + 7r). The Josephson current exhibits a jump when a pair of 
Andreev bound states aligns with the Fermi energy. The condition for the current jump is given. In 
particularly, we find that the position of the current jump and the position of the maximum value 
of the critical current I c are identical. Due to the interference between the two paths, the critical 
current I c versus the QD's level Ed shows a typical Fano shape, which is similar to the Fano effect 
in the corresponding normal device. But they also show some differences. For example, the critical 
current never reaches zero for any parameters, while the current in the normal device can reach zero 
at the destruction point. 
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PACS numbers: 73.63.Kv, 73.23.-b, 74.45.+C 



I. INTRODUCTION 

Mesoscopic electron transport through an Aharonov- 
Bohm (AB) interferometer has attracted considerable 
attention recently because of its applications in nano- 
technology. It reveals information about intrinsic quan- 
tum states by detecting interference of electrons in dif- 
ferent paths. Interference between a continuum energy 
spectrum and a discrete energy state gives transmission 
probability T(E) asymmetric line shape of typical Fano 
resonance. 1 Proposed first by Fano, it has been broadly 
studied and observed in recent experiments.—^ By set- 
ting a quantum dot (QD) in one arm of the AB interfer- 
ometer, the Fano resonance and Kondo-Fano resonance 
are founds Interference between direct transmission and 
QD makes the transport phase through QD observable^ 
In addition, the extra phase <I> due to cither the mag- 
netic flux or the spin-orbital interaction can modulate the 
interference^ The transmission probability T(E) shows 
periodical function of the extra phase. In such an AB- 
Fano system, the T(E) versus QD level shows a Brcit- 
Wigner resonance when the phase $ is tt/2 and a typical 
Fano resonance when the phase is or ir. The T(E) also 
shows a Breit-Wigner resonance when the direct tunnel- 
ing is broken and a Fano type when direct tunneling is in- 
creased. Both the phase parameter and direct tunneling 
can be included in a Fano parameter q with the transmis- 
sion probability T(E) cx ^j-, in which e = E ~ r £d , and 
E is the energy of the incident electron, Ed is the QD's 
level and T is the coupling between QD and the lead. The 
Fano parameter q, generally a complex number, decides 
the line shape of the resonance. 

Previous work on the AB-Fano interferometer focused 
mainly on the normal device, in which the two exter- 
nal leads are normal. By attaching the interferometer to 
two superconductor leads instead of two normal leads, 
the Cooper pair transport and Andreev tunneling occur, 



then Josephson current emerges even at the zero bias. 
The purpose of this paper is to study how the Josephson 
current is affected by the interference of the two paths, 
and whether the Josephson current also has the Fano 
characters as in normal systems. Recently, the Joseph- 
son current through a mesoscopic system has been ex- 
tensively investigated because of scientific interest and 
possible applications. The Josephson current through a 
clean thin superconductor-normal-superconductor (S-N- 
S) junction has a discontinuous jump at superconduc- 
tor phase difference 9 = tt under proper conditions. 8 
This current jump arises from the discontinuity of su- 
percurrent contribution of Andreev bound-states driven 
by phase difference 9 of coupled leads.— In this situa- 
tion, the current jump can be broken by impurities, 10 
the finite temperature, the normal lead attachment^ or 
the electron-electron interactional- Another interesting 
characteristic of the Josephson current is the 7r junction, 
where the sign of Josephson current can be reversed from 
I c sin(0) to I c sin(# + ir)r^ The tt junction is broadly re- 
searched in S-ferromagnetic-S (S-F-S) junctions , 14 ! 15 S-F 
ringi£ S-QD-S junctions ,±2^ S-AB-S junctions,-^^ anc j 
so on. 

In this paper, we will investigate the Josephson cur- 
rent / through an AB-Fano interferometer consisting of 
a QD in one arm and a magnetic flux $ through the ring. 
By using the non-equilibrium Green's function method, 
the Josephson current expression is obtained. Due to the 
electron-hole symmetry, the Josephson current has the 
property 7(£<2, = I{—£d,& + 7r). Without the mag- 
netic flux, the current-superconducting phase (1-9) rela- 
tion usually shows a sinusoidal line shape. But when the 
magnetic flux exists, the curve of 1-9 is quite complex 
and can be modulated to tt junction in appropriate pa- 
rameters. The system has two pairs of Andreev bound 
states due to the direct arm and the QD. When one pair 
of Andreev bound states is in a line with the Fermi energy 
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Ep = 0, the Josephson current jumps. The conditions 
for this current jump are given. In particular, we find 
that the position of the current jump is identical to the 
position of the maximum value of the critical current I c . 
The critical current I c versus the QD level E d shows a 
typical Fano shape due to the interference between the 
two paths. The positions for the constructive and de- 
structive interferences are same with the corresponding 
normal device. But the critical current can not reach 
zero at the destructing position, which is different from 
the normal device in which the current can be zero at 
the destructing position when the magnetic flux $ = rnr 
(integer n). In addition, the critical current is a periodic 
function of the magnetic flux with the period ir at the 
QD level Ed = and period 2ir while Ed ^ 0. 

The rest of the paper is organized as follows. In Sec. 
II, the Hamiltonian is present and the Josephson cur- 
rent expression is derived. Main numerical results are 
given in Sec. Ill, in which we investigate the Josephson 
current-superconducting phase relation, the condition of 
the jump of the Josephson current, and the characters of 
the critical current. A brief conclusion is given in Sec. 
IV. 



Hamiltonians are transformed to the following forms: 

Ha — 22 EkCl a a Ckcr,a 
ka 

k 

H T = ^{t a e^Cl aa d a + t a e^ SL dlCka, a ) 

k(T,Oi 

k,k' ,er 

Hdot = ^Edd^da. (1) 



with <p = 2$ - 8 L + R . 

The current I through AB-Fano interferometer can be 
calculated from the evolution of the electron number op- 
erator N L = J2 k * C Ll C ^l in the left leads^i 

I = ~e(N) (2) 
= yRe J ^[t L e^G< LA1 (E)+te-^G< LA1 (E)] 

Here the Nambu representation has been used. The 
Green's function G dL {E) and G RL (E) are the Fourier 
transformation of G dL (t — t') and G RL (t — t'), which are 
defined as: 



II. MODEL AND FORMULATION 



The system we considered is the AB-Fano interferom- 
eter consisting of a QD and a reference arm connected 
to two BCS superconductor leads. The Hamiltonian of 
the superconductor leads is H a = J2ka e kC\ a a Cka,a + 
£ fc (A Q C fei , Q C_ feT , tt + A* a Cl kha Cl , J, where a = 
L, R represent the left and right lead, and A Q = 
Ae l8 ° is the complex superconducting order parame- 
ter, with the superconductor gap A and the supercon- 
ductor phase a . Coupling Hamiltonian between leads 
and QD is, H T = ^Z^JtaP^JL,, + t a d\C ka ^ a ) + 

T, k ,k'A tRLC l-y,L C k'<yM + tLRCl R C ka>L ). Parameter 
t a is the coupling coefficient of the QD to a-th lead, 



and t LR = t* RL = te l * = 



te *° 



is the direct cou- 
pling between the two leads. With the consideration of 
the magnetic flux <f> through AB interferometer, a phase 
<f> = 2ir4>/4>Q is added to the coupling coefficient t, with 
4>o = e/h the flux quantum. We adopt the single level 
QD with neglecting the intra-dot Coulomb repulsion; its 
Hamiltonian is Hdot — 53<t £ dd\d a . The total system 
Hamiltonian is H = Hl + H R + Hdot + Ht- 

For the sake of calculation, we take an unitary transfor- 
mation with U = exp{J2 kt7 , a % \ G ka, a C ka,a\? 2 and the 



W <^ Ti (tW)> (c^ L (t')d\(t)) 



v ( (cl, L (t')c k ' hR (t)) (ct r , L (t')cl k , UR (t)) 

^ \ (c- kltL (t')c k , hR (t)) (c- kL L(nci k , LR (t)) 



fr.fr 



Notice here that all of the Green's functions are functions 
of time difference t — t' , because that we investigate the 
dc Josephson current in the zero bias case. 

Also, because of the zero bias case, the system 
is in equilibrium and the fluctuation-dissipation theo- 
rem holds, we now have G< L (E) = -f(E)(G r dL {E) - 
G a dL {E)) and G< L {E) = -f(E)(G RL (E) - G RL (E)), 
where f(E) = l/{exp[(J3 - E F )/k B T] + 1} is the Fermi- 
Dirac distribution. The current expression becomes 



1 - t / t ! ^ G <- 



G dL) f L 



i*(Gk-G^)M3) 



In the following, we need to solve the retarded 
Green's function G r dL (E), G r RL (E). In the Nambu 
representation, the retarded Green's function of iso- 
lated superconducting leads and QD are respectively^ 

Sa (E) = and 8dd (E) = 

l/(E-£ d + i V ) 



l/(E + e d + it}) 



, where p{E) is the 
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normal density of states, r\ is an infinitesimal real num- 
ber, (3 (E) = pA/E, and 0(E) = 
\E\ < A and 0(E) = i^j^fE 2 
A. Tunneling coefficients in the Nambu representa- 
tion can also be expressed in 2 x 2 matrix with: t a = 
e ie a /2 \ / elv /2 o 



E/^fA 2 - E 2 while 

^A 2 while \ E \ > 



t, 







-ie a /2 



and t LR = t 







In the following calculation, we take the symmetric barri- 
ers with tL = t R for convenience. By using Dyson's equa- 
tion, the Green's function g r of the system decoupling 
with the QD (i.e. tL =t R = 0) can be deduced as: g£ L = 

(gr 1 - t M gkt^)-\ g RR = (g^- 1 - t* LRg it LR )-\ 

and g r LR = g r L t LR i r RR = g r LL t LR g R . Then the re- 
tarded Green's function of the whole AB-Fano interfer- 
ometer device are solved as: G r dL — G^ d (t£g£ +t* R <g RL ), 
G r RL = ¥ RL + &R + ¥ RL t L )G r dd {tl¥ L + t* R g RL ), and 



dd 



g. 



dd 



where the retarded self energy S r 



G 

is: ST = t* g£t L + t* R g R t R + t* L gl R t R + t* R g RL t L . 

After solving the retarded Green's functions G r , the 
expression of the Josephson current / through AB-Fano 
interferometer can be reduced as: 

-prdE 

h / 2tt 3 y ' 1 V 



VP 



(QiiA u + Q 12 Ai 



(4) 



T^p 2 e~ i — 



Qn + A\ 2 Q 22 + 2A U A 12 Q 12 }} 



where V = 1 + x 2 + 2x(E 2 - A 2 cos ip)/(E 2 - A 2 ) and 



arm and QD. The Josephson current in the equation (4) 
can be split into two parts, the continuous part I con , 
contributed from continuous spectrum while the energy 
E outside the superconducting gap A, and the discrete 
part Idis contributed by Andreev bound states while E 
within the gap. In numerical calculation I con is obtained 
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FIG. 1: Current / vs. superconductor phase 6 for various 
parameters: (a) $ = 0, x — 0, and F = 0.1 for different QD's 
level Ed', (b) $ = and r = for different x; (c) and (d) 
x = T = 0.5, 4> = (c) and ty (b) for different Ed- 



E = \EV + rp(l + x)} 2 - [edV-TyGix- P 2 )cos— 
- T^pl cos B —^\ 2 - r 2 /3 2 [cos - + xco^- + ^)] 2 (5) 

Here x = t 2 Tr 2 p 2 dictates the tunneling through the di- 
rect arm and T = 2irpt 2 a is the coupling strength of QD 
to leads. In the wide-band approximation, x and T are 
independent with the energy E. The factors A and Q in 
equation (4) are: 



(l + x) + ^(x + P 2 )e 

A 



An 

Al2 

A 21 = #e-*V 2 (l + 
A 22 



R2 







E e i9 ' 2 {l + xe*P) + 2isin§ § /3^i 



(l + .T) + ^(a; + /3 2 )e- 



— 2i sin 



Po e 



by integral in equation (4) and Idis is approached by 
dealing with the delta functions due to the infinitesi- 
mal imaginary part irj. The discrete part is usually 
much larger than the continuous part. The factor 3 
has two pairs of poles at E^ 2 (E^ 2 = —E^ 2 ) within 
the gap, which count for Andreev bound states. These 
Andreev bound states arise from the hybridization of 
bound states of QD Eq D and the direct arm Ef. Here 

Eq = ±AyJl — pj^p sin 2 (y>/2) are the Andreev bound 

states in the direct arm^ and Eq are also the poles of 
the factor T>. Eq D {Eq D = —Eq D ) are the Andreev 
bound states in the QD for the S-QD-S device. 



and 



On = EV + e d V - TpRe(A 22 ) 

0i 2 = TpRe(A 12 ) 

Q 21 = TpRe(A 12 ) 

Q 22 =EV- s d V + TpRe(A n ) 



The first term in equation (4) describes the direct 
tunneling contribution, the second and third terms are 
transport through QD and interference between direct 



III. NUMERICAL RESULTS 



In this section, we present our numerical investigations 
on the relation of the Josephson current versus the su- 
perconducting phase, the condition of the jump of the 
current, the Fano resonant characters of critical current, 
and the dependence of critical current on the magnetic 
flux. 
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Josephson current-superconducting phase 
relations 



viously wider than the positive-current region. However, 



We first discuss the current-phase {1-9) relation of AB- 
Fano interferometer. Fig.l shows Josephson current I as 
a function of the superconducting phase difference 9 of 
the left and right leads. When x = 0, the device reduces 
into a QD coupled to two superconductor leads. The 1-9 
relation is sinuous line shape when the level Ed is far off 
the Fermi level Ep = 0, and the current shows a discon- 
tinuous jump at 9 = tt when Ed — (shown Fig. la), in 
which the Andreev bound states Eq D = Eq D = 0£L2££^ 
On the other hand, while r = 0, the system reduces into 
an S-FS device. The 1-9 relation is sinuous line shape 
when x is much smaller or larger than 1 , and the current 
/ has a discontinuous jump while x = 1 (see Fig. lb), 
in which the Andreev bound states Eq = Eq = . 10 i 26 
When both x and T are non-zero (in other words, the 
two pathes are opened), the AB-Fano interferometer is 
formed and the transport can be adjusted by magnetic 
flux $ through it. Two pairs of Andreev bound states 
and Eq D , which belong to the direct arm and QD, cause 
the hybridization to form new Andreev bound states 
E^ 2 i which enable the interference construction or de- 
struction of the Josephson current. Then the current- 
phase 1-9 relation is usually not a sinusoidal- like function 
(except for the special magnetic flux values <& = and 7r) . 
While $ = and 7r, the 1-9 relation is still a sinuous line 
shape, and the current / is zero at 9 = and 7r, as shown 
in Fig.lc and Id. In some specific parameters the dis- 
continuous jump of the current can still occur, which we 
will detail in the next sub-section. Here, we notice that 
the current has the relation: I(£d, = I(—Ed, $ + 71"). In 
other words, while the level changes from £d to — Sd and 
the magnetic flux from <E> to <& + tt, the current / does 
not vary regardless of any other parameters. The rela- 
tion of I(ed, $) = I(—Ed, $ + 7r) comes from the electron- 
hole symmetry, i.e. taking the transform (d a ,dka-,a) to 
(JJ., dj^ Q ) an d simultaneously setting the parameters 
(ed, 3>) to (—Ed, $ + 7r), the Hamiltonian H is invariable. 

When the magnetic flux $ is not equal to or tt, 
the current-phase relation is usually not a sinusoidal- 
like function, and the current / has non-zero values at 
9 = 0. In some special parameters, the current / is nega- 
tive while 9 £ [0, tt], which is a tt junction. For example, 
by proper selection of parameters, Ed = 0.5, T = 0.45, 
<3> = 0.67T, and x = 0.4 ~ 0.9, the current / is negative 
when the phase 9 £ [0, tt] as shown in Fig. 2a. But it is 
not a strict tt junction, and the current is not positive in 
all regions 9 £ [tt, 2tt]£2- The realization of quasi- tt junc- 
tion is because of the introduction of the magnetic flux 
phase 4> which changes the interference of two pathes. 
Another example, as shown in Fig. 2b with the magnetic 
flux $ = tt/2, the current is negative (positive) in most 
parts of region 9 E [0,tt] ([tt, 2tt\). Here we notice that 
the 6 region for the negative current / is not equal to 
that of the positive current. In all curves in Fig. 2a and 
some curves in Fig. 2b, the negative-current region is ob- 
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FIG. 2: Current / vs. the phase 9 with the parameters in (a) 
£d=-0.5, r=0.45, <!>=0.67r, and x from 0.4 to 0.9 with space 0.1 
in arrow direction, (b) £=0.5, T = 0.1, <&=0, and Ed = —0.1 
(solid curve), —0.05 (dashed curve), (dotted curve), 0.05 
(dash-dotted curve), and 0.1 (dash-dot-dot curve). 

no matter the parameters, the positive-current region al- 
ways exists. In other words, it is not possible for the 
current to be negative in the entire region 9 £ [0, 2tt]. In 
addition, while Ed — 0, a jump emerges in the curve of I- 
9 at 9 = tt (see Fig. 2b). We will study this discontinuous 
jump in detail in the next sub-section. 



B. the condition of the jump of the current 

The current I con from continuous spectrum with \E\ > 
A is always continuous. The discontinuous current arises 
from the part Idi S which is from the Andreev bound 
states. When one of the two pairs of Andreev bound 
states E^ 2 just aligns with the Fermi level Ep = (i.e. 
E+ = E^ = 0), an abrupt jump occurs in the current 
Idis so that the current I = I con + Idis- The condition of 
the jump of the current / is thus = 0) — 0. With 
the help of Eq.(5), the condition S(_B = 0) = can be 
reduced to: 



£ d T> ~ Txy/x cos® ~ Ty/x cos(9 - = 
r 2 [cosf +xcos(2$-6y2)] = 



(6) 



When the above two equations are tenable, the current 
will jump. For example, (i) while T = 0, the condition in 
Eq.(6) reduces into x = 1 and 9 = tt. This is consistent 
with the jump in Fig. lb. (ii). While x — 0, the condition 
in Eq.(6) reduces into Ed = and 9 — tt, which agrees 
with the jump in Fig. la. In general, when both x and T 
are non-zero with a magnetic flux through the AB-Fano 
interferometer, the condition in Eq.(6) can be rewritten 
as: 

r l + a ; cos2$ 1 

9 = -2 arctan . ^ (7) 

1 xsh2$ J w 
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FIG. 3: Upper panel shows the jump position in the pa- 
rameter space of (x, £d/r). The curves from left to right 
are for $ = 0.17T, 0.3-7T, 0.7n, and 0.9-7T. Lower panel shows 
the current / vs. the QD's level Ed with x — 0.2 and 
9 = — 2arctan[(l + a;cos2$)/(a;sin2<I>)]. The magnetic flux $ 
for different kind of curves are same with the upper panel. 



the same as the above superconducting device. From 
the Eq.(9), we found that the interference construction 
can be realized at 




FIG. 4: Critical current I c vs. Ed with parameters in (a) 
x=0.l and $=0; (b) £=1.0 and $=0; (c) and (d) F = 0.5 and 
$ = 0; (e) T = 0.5 and x = 0.2; (f) T = 0.5 and x = 0.8. 



and 



Ed = — 



r-^/xcos^ 

1 - X 



(8) 



Specifically, for $ = 7r/2, we have 9 = ir and = 0, 
which is consistent with the jump position in Fig. 2b. 
Fig. 3a shows the jump-occurrence region in the param- 
eter space of (x, ^). When the parameters are just at 
the curves of Fig. 3a, a current jump will occur. Fig. 3b 
shows the current / versus the level e^/r at x = 0.2 and 
r = 0.5. It clearly shows that the jump occurs at the 
corresponding intersections in Fig. 3a. 



C. Fano resonant characters of critical current 

In this and the next sub-section, we will focus on the 
critical current, which is experimentally accessible. In 
fact, the critical current in the superconducting AB-Fano 
interferometer behaves similarly to the current in the nor- 
mal AB-Fano interferometer under the small bias. So 
here, we simply recall the results of the normal AB-Fano 
device, which consists of an AB ring attached to two nor- 
mal leads with a QD in one of its arms. The transmission 
probability T of the normal AB-Fano device is^£ 



T(E) 



Ax 4r(l - x)y/xcos& 
(1 + xf + {l + xf 
r[(l + a;) 2 -4a;(l + cos 2 $)] 



ReG r (i;) 



(l + x) 3 



ImG r (£), (9) 



where G r -\E) =E-e d + T^cos $/(l +x) - iT/(l + x) 
and the meaning of the parameters (sd, T, <&, and x) is 



Y^fx cos $ 



and interference destruction is at 



(1 + .TCOs2$)r 
2 % /i(l + x)cos<I>' 



(10) 



(11) 



Here the interference construction position is the same 
as Eq.(8) of the position of the jump of the supercon- 
ducting current in the superconducting device. The I- 
Ed relation of the normal device shows a symmetric line 
shape when x — and a typical Fano resonance when 
at finite x and zero magnetic flux $. The shape of the 
curve of I-Sd can be affected by the magnetic flux $ and 
is symmetric at $ = ir/2& 

Next, we study the critical current I c in the super- 
conducting AB-Fano interferometer. Here we select the 
maximum of current in a period 27r of superconducting 
phase difference 9 as the critical current. Due to the in- 
terference between two pathes, the construction or de- 
struction transport occurs, and the critical current I c 
versus the QD's level Ed usually exhibits a Fano char- 
acters. To show the details of Fano characters, we plot 
I c as a function of Ed at zero magnetic flux (<I> = 0) for 
small x values and large x values in Fig. 4(a) and (b), 
respectively. For small x values (x = 0.1, for instance), 
I c shows a Fano resonance near Ed = with the peak 
at Ed = — T^/x/(l — x) and the valley at approximately 
Ed = T/(2y/x). The positions of the peak and valley are 
the same with the normal device [see Eqs.(10) and (11)], 
though in the normal device the normal current is un- 
der the small bias, while in the present superconducting 
device the critical current I c is at zero bias. With the 
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enhancement of T, the line shape tends to behave more 
symmetric characters as the Breit-Wigner resonance and 
the magnitude of the critical current is enhanced at neg- 
ative Ed side but reduced at positive Ed side. For large x 
values (e.g. x = 1.0, as shown in Fig. 4b), the destructive 
interference plays the core role, and the curves of I c -£d 
are almost symmetric for the small T and show an obvi- 
ous Fano valley at Ed = T/2. The Fano peak is pushed 
to infinity at x — 1.0 (as shown in Eq.(10)), so it is in- 
visible in Fig. 4(b). With T increasing, the transmission 
probability through the QD grows and the destruction of 
two paths occurs, causing the valley to deepen and the 
line shape to become more asymmetric. 



In this paragraph, we investigate the effect of the direct 
path x on the critical current I c . In Figs. 4(c) and 4(d), 
we show the critical current I c versus the QD's level Ed 
at different x for a fixed T — 0.5 and $ = 0. While x = 0, 
the direct path is closed, I c is completely symmetric with 
the peak at Ed = and the valley at Ed = oo, which are 
the same with Eqs.(10) and (11). When x increases, the 
peak position is moved off from the Fermi level Ep = 
but the valley position gradually from infinity to 0, and 
the curve of I c -£d shows a Fano resonance. At the middle 
x, the Fano resonance is the most prominent. While x 
tends to 1, the direct path is completely open and the 
peak position tends to infinity, leading the I c -£d curve 
into a symmetric valley. 



In the above I c -£d relation discussion, the magnetic 
flux <f> is fixed at zero. In the normal AB-Fano device, 
the line shape of the current I versus Ed is strongly mod- 
ified by the magnetic flux $. So, in the following, we 
investigate how the I c -£d curve in the superconducting 
AB-Fano device is affected by Figs. 4(e) and (f) show 
I c -£d with r > x and T < x, respectively. These two 
cases when T > x and T < x represent situations where 
the QD path or direct arm path dominates the trans- 
port. While r > x, the I c -£d curve is strongly affected 
by the magnetic flux <f>. When $ increases from 0, the 
Fano shape first grows more notable, then the peak is 
greatly reduced while $ near n/2, and at last the peak is 
increased and the Fano shape is recovered again at $ = n 
(see Fig.4e). On the other hand, while T < x, the I c -£d 
curve is only slightly affected by $ (see Fig.4f). We add 
the following three points: (i) The critical current I c has 
the relation: I c (sd,^) — I c (—£d,& + tt), because of the 
electron-hole symmetry and /(£<;, <f>) = I(— £d,$ + 7r)- 
Due to the relation I c (£d, $) = I c (—£d, $ + n), the criti- 
cal current I c is a periodic function of $ with period n at 
Ed = 0. (ii) At $ = 7r/2, the I c -£d curve is still not sym- 
metric, which is different from the current in the normal 
device, (iii) Though the valley in the normal device can 
reach zero in some parameter regions, at no parameters 
does the valley value of the critical current I c reach zero. 




0.0 0.5 1.0 1.5 2.0 



0(71) 



FIG. 5: Upper panel: critical current 7 C vs $ for different x 
with Ed — and Y = 0.5. Lower panel: critical current I c vs. 
<£> for different e d with x = 0.2 and V = 0.5. 



D. Critical current-magnetic flux relations 

Finally, we investigate the relation of the critical cur- 
rent I c with the magnetic flux phase $. Fig. 5(a) shows 
I c versus $ for the different x values at Ed = 0. Several 
characteristics can be noticed: i) I c is a periodic func- 
tion of $ (or the magnetic flux (j>) with the period 2ir (or 
e/h) while Ed ^ 0, and period ir (or e/2h) at Ed = 0. ii) 
For a small x value, I c is almost a constant, because the 
direct path is almost closed. For a large x value, the os- 
cillation of I c -<f> is strong. The oscillation shape greatly 
departs from sin <!> or cos $ shape because the higher or- 
der tunneling processes are numerous at Ed = 0. iii) 
The critical current I c shows a peak at $ = (2n + l)n/2 
with the integer n, which is consistent with equations 
(7) and (10). Fig. 5(b) shows I c versus $ for the QD 
level Ed at x = 0.2. When Ed is far away from zero, 
I c shows a standard sin<J>/cos$ behavior because the 
higher order tunneling processes are weak at the small 
x and \£d\ S> 0. But I c for positive Ed bears a phase 
lapse of it with that of —Ed- This phase lapse is from the 
phase Qqd of the transmission coefficient of a QD, with 
Qqd = tt/2 when > and — tt/2 when Ed <§C 0. As 
Ed approaches zero, the oscillation of I c -<& is enhanced 
and departs from the sin<I>/ cos$ behavior, with peaks 
emerging at about <f> = 7r/2 and 37r/2. 



IV. CONCLUSION 

In conclusion, the Joscphson current through an 
Aharonov-Bohm interferometer consisting of a quantum 
dot and a direct arm with magnetic flux through the ring 
has been investigated. An equation for the occurrence 
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of supercurrent discontinuity is given. In particular, we 
found that the position of the supercurrent discontinuity, 
the position of the peak of critical current, and the con- 
structive interference of the current in a corresponding 
normal device are the same. By adjusting the device's pa- 
rameters, such as the magnetic flux phase and direct arm 
coupling, the Josephson junction can vary from a normal 
junction to a 7r-junction. Fano characters of the critical 
current are similar to those of the current in the normal 
device in the small bias, but there are also some differ- 
ences. For example, the critical current cannot reach zero 
with any parameters, while the current in the normal de- 
vice can reach zero in the destruction position. Also, the 
critical current with variation of magnetic flux shows a 



period of e/2/i when QD level aligns to the Fermi level, 
or a quantum magnetic flux of e/h when the QD level 
apart from the Fermi level. 
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